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One dimensional free-fermions and hard-core bosons are often considered to be equivalent. Indeed, 
when restricted to nearest-neighbor hopping on a chain the particles cannot exchange themselves, 
and therefore hardly experience their own statistics. Apart from the off-diagonal correlations which 
depends on the so-called Jordan- Wigner string, real-space observables are similar for free-fermions 
and hard-core bosons on a chain. Interestingly, by coupling only two chains, thus forming a two- 
leg ladder, particle exchange becomes allowed, and leads to a totally different physics between 
free-fermions and hard-core bosons. Using a combination of analytical (strong coupling, field the- 
ory, renormalization group) and numerical (quantum Monte Carlo, density-matrix renormalization 
group) approaches, we study the apparently simple but non-trivial model of hard-core bosons hop- 
ping in a two-leg ladder geometry. At half-filling, while a band insulator appears for fermions at 
large interchain hopping t± > 2t only, a Mott gap opens up for bosons as soon as tx 7^ through a 
Kosterlitz-Thouless transition. Away from half-filling, the situation is even more interesting since a 
gapless Luttinger liquid mode emerges in the symmetric sector with a non-trivial filling-dependent 
Luttinger parameter 1/2 < Ka < 1. Consequences for experiments in cold atoms, spin ladders in a 
magnetic field, as well as disorder effects are discussed. In particular, a quantum phase transition 
is expected at finite disorder strength between a ID superfluid and an insulating Bose glass phase. 

PACS numbers: 05.30. Jp,75. 10. Jm,05.30.Rt 



I. INTRODUCTION 

Low-dimensional bosonic systems have attracted an 
increasing interest during the last few decades. The 
competition between various ingredients such as interac- 
tions, quantum fluctuations, geometrical constraints (low 
dimensionality/coordinance, frustration), possibly disor- 
der, may lead to a large variety of interesting states of 
matter, as for instance the enigmatic supersolid phase 
of *He under pressure.^ Perhaps the most extreme situ- 
ation arises for one-dimensional (ID) systems where in- 
teractions are known to induce dramatic effects. While 
ID systems have long been considered as purely abstract 
objects reserved for theoretical studies, they have now 
become experimentally accessible. In particular, due to 
a very intense activity, several systems are now available 
in the lab to achieve ID bosonic physics in various con- 
texts: Josephson junction arrays, ^"^^ superfluid ^He in 
porous media, ^^"^^ cold atoms, ^'^^"^^ quantum antifer- 
romagnets.^^""^^ Furthermore, several theoretical studies 
have recently shown that ID physics is still (80 years af- 
ter Bethe's seminal work)"^^ a very active playground for 
exploring novel and exotic phases of matter. •^^■•^''"''^ 

As in most condensed matter problems, only a few 
competing terms dictate the physics of ID bosonic sys- 
tems. For instance, when kinetic processes dominate over 
(repulsive) interactions, a delocalized superfluid (SF) 
phase is favored. Note that despite the absence of a true 
Bose-condensed phase in ID systems,''^ superfluidity is 
still possible at zero temperature since it only requires 
off-diagonal quasi long-range order, the resulting phase 
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FIG. 1: (Color online) Schematic picture for the two-leg lad- 
der system. 



being a Luttinger liquid (LL). Conversely, when strong 
repulsion dominates, a Mott insulator (MI) is expected, 
either breaking or not lattice symmetries, depending on 
the commensurability of the particle filling. 

An interesting case in ID is the hard-core boson limit 
(also known as the Tonk-Girardeau g achieved in 

cold atom experiments)^^ which, on a lattice, corre- 
sponds to infinite on-site repulsion. In such a limit, and 
in the absence of density-density interactions between 
neighboring sites, an exact mapping between fermions 
and hard-core bosons'^^ implies that bosons behave "al- 
most" like free-fermions, at least regarding their real 
space properties, where a "real space Pauli principle" 
holds. However, when the lattice geometry allows for 
particle exchange, the bosonic statistics may play a ma- 
jor role, leading to qualitatively different phase diagrams. 
It is quite clear in two dimensions, ^^■'^^ but for quasi-lD 
geometries like ladders one may wonder whether a min- 
imal model made of two coupled chains would allow the 
bosonic statistics to emerge and lead to qualitatively dif- 
ferent physics as compared to free-fermions. 

In this paper, we aim at exploring one of the sim- 
plest quasi-lD model where free-fermions and hard-core 
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bosons display qualitatively distinct phase diagrams. 
The model, defined on a two-leg ladder (depicted in 
Fig. 1), is governed by the following Hamiltonian 

H ^ -t NjVj+i -A* "^^'3 

-*-LE[^L^2,,+h.c.], (1.1) 
j 

where the sites along each leg £=1,2 are labeled by 
j = 1 . . . L, L being the total length of each chain. The 
operator b\ ■ creates a particle on site and t, t± 

are the longitudinal and transverse hopping integrals re- 
spectively. Tiij = bl jbgj is the onsite density opera- 
tor and /i is the chemical potential. The filling of the 
ladder will be denoted by p = N/ {2L) in the following, 
with N the total number of bosons. The main phases of 
(1.1), discussed in this work, are summarized on Fig 2 
in the cases of free-fermions (left) and hard-core bosons 
(right). For non- interacting fermions, the phase diagram 
is trivially obtained by filling two bands: the system is 
a simple metallic state at all fillings < p < 1, except 
when the interchain hopping t± > 2t where a band in- 
sulator appears at half-filling p = 1/2. The situation is 
very different with hard-core bosons for which there is 
no Pauli principle when filling the states in momentum 
space. Indeed, as demonstrated below, the interchain 
hopping t± is shown to be a relevant perturbation at 
half-filling, opening up a charge gap ~ exp{~at/t±) 
via a Kosterlitz-Thouless mechanism. Therefore, in con- 
trast with frec-fermions, the bosonic insulating state at 
half-filling, called rung-Mott insulator, extends down to 
= 0. For incommensurate fillings, the bosonic sys- 
tem is a single mode LL (called LL^-'^^ by contrast with 
LL^^) for the two-band or two-mode Luttinger liquid of 
free-fermions, see Fig. 2) with a finite superfluid (SF) 
fraction and, most interestingly, a strongly renormalized 
Luttinger parameter K{p) which varies continuously be- 
tween 1/2 and 1. 

In the rest of the paper, using a combination of ana- 
lytical and numerical techniques, we discuss in details all 
these non-trivial aspects. In section II, after recalling the 
similarities between free-fermions and hard-core bosons 
on a chain, we present some well-known exact results for 
free-fermions on a ladder. In section III, the two coupled 
chains problem is studied for hard-core bosons in the two 
analytically accessible limits: strong and weak interchain 
couplings, using a perturbative approach, bosonization 
and renormalization group (RG) techniques. We also 
compare the physics of this model with the well known 
case of quantum spin ladders. We then use exact numeri- 
cal techniques, quantum Monte Carlo and density-matrix 
renormalization group, and confront them with analyt- 
ical predictions in order to provide quantitative results 
for several quantities of importance: the bosonic densi- 
ties (normal and SF), the charge gap, correlation func- 
tions, correlation lengths, and Luttinger parameters. In 
particular, section IV focuses on the half-filled insulating 
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FIG. 2: (Color online) Phase diagram of the two-leg ladder 
model Eq. (1.1) for free-fermions (left) and hard-core bosons 
(right). The metallic phase for fermions comprises either one 
LL''^' or two bands LL'^-*, whereas the bosonic SF is a sin- 
gle mode (symmetric, see the text) Luttinger liquid LL'-^^ 
For hard-core bosons, the rung-Mott insulating phase extends 
down to = 0. 

state while the LL behavior at incommensurate fillings 
is studied in great details in section V where the global 
phase diagram is shown. Finally, in section VI we dis- 
cuss two important issues related to experiments in cold 
atoms and spin ladder materials as well as the effect of 
disorder on this bosonic ladder. Section VII concludes 
the paper. 



II. FREE-FERMIONS VS HARD-CORE BOSONS 

A. Quantum statistics in ID systems 

The concept of quantum statistics stems from the 
indistinguishability of identical particles in many-body 
quantum systems. If 5'(xi,X2) is the wave- function of 
a system made of two indistinguishable particles, the 
probability of finding two particles at position xi and 
X2 must be invariant under the exchange or particles, 
that is, |^'(xi,X2)P = |5'(x2,xi)p. Therefore ^'(xi,X2) 
and \E'(x2,Xi) are equal, up to a phase factor e*". In 
three dimensions, it can be argued that a further ex- 
change of particles is actually equivalent to the identity 
transformation, imposing that a be an integer multiple 
of TT. Two situations arise: the wave-function is either 
symmetric or antisymmetric under the exchange of two 
particles, thus defining bosonic and fermionic statistics, 
respectively. Note that the present restriction does not 
hold in lower dimensions where so-called fractional statis- 
tics are expected to exist. However, even though we will 
only be interested in ID and quasi ID systems we will 
focus on bosons and fermions only. This real-space ap- 
proach of quantum statistics extends to any Hilbert space 
in that a given many-body state - a Fock state - will 
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either be totally symmetric or totally antisymmetric un- 
der the exchange of any two particles. As a consequence, 
fermions obey Pauli's exclusion principle, while several 
bosons can indeed occupy the same quantum state - ac- 
tually a macroscopic number in the phenomenon of Bose- 
Einstein condensation. There is one situation though, 
where differences between bosons and fermions are less 
drastic. Indeed, a ID gas of impenetrable bosons ex- 
hibits similarities with a ID gas of spinless fermions. 
The hard-core constraint mimics the exclusion principle, 
while the ID character ensures that the particles cannot 
move around each other. Wave-functions are identical up 
to a symmetrization factor enforcing the quantum statis- 
tics. Therefore, quantities depending only on the mod- 
ulus of the wave-function, such as density correlations, 
are identical, whereas off-diagonal quantities, such as the 
momentum distribution, are affected by statistics. These 
various quantities - density and momentum distribution 
- are readily encoded in the single-particle density-matrix 
{tp^ {x)ip{x')), where ip^{x) and ipi^) ^re creation and an- 
nihilation operators. 
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FIG. 3: (Color online) Schematic representation of the mo- 
mentum distribution for ID free- fermions (left) and hard-core 
bosons (right). 

in the continuum case, there is a one to one correspon- 
dence between the eigenstates of the Fermi Hamiltonian 
and the eigenstates of the Bose Hamiltonian. 



When fermions and hard-core bosons look alike 



2. When fermions and hard-core bosons look different 



In the continuum, if f{xi, . . . ^ xn) is a solution of 
the Schrodinger equation for a ID gas of N fermions, 
satisfying the constraint ^'i?(xi, . . . , xat) = when- 
ever Xi — Xj for all i and j, then ^^(a;!, . . . , a;jv) = 
A{xi, . . . jXn)"^ f{xi, . . . ,xn) is also a solution of the 
Schrodinger equation for a ID gas of hard-core bosons, 
with A = Yiiyj^^S^i^i ~ antisymmetric func- 

tion which can only be equal to ±1.'^^ It appears clearly 
that f{xi, ■ ■ ■ ,Xff)\^ — |\E's(a;i,...,a;7v)P, so that any 
quantity depending only on the modulus of the wave 
function is identical for fermions and hard-core bosons. 
As noticed by Girardeau, "^'"'^'^ this one to one correspon- 
dence between fermionic and bosonic wave-functions is 
only possible in one dimension. The hard-core constraint 
divides the parameter space in iV! disjoint regions, and A 
changes sign discontinuously at the boundary of each re- 
gions. In two and three dimensions, the parameter space 
remains connected and a function equivalent to A cannot 
be defined. On a ID lattice the correspondence is ensured 
through the Jordan- Wigner transformation.^^ Hard-core 
bosons on a lattice have mixed commutation relations: 
[h,b]] = = [bib]] = for z ^ J, {6„6l} = 1 and 

{bi)"^ = {b})'^ — 0. The Jordan- Wigner transformation 
maps this problem of hard-core bosons onto a problem 
of spinless fermions c; through: 



i-i 

ci = exp I i-K b]bj 



(2.1) 



The Jordan- Wigner string, although making the transfor- 
mation non-local, ensures that c; and c[ satisfy anticom- 
mutation relation and are indeed fermionic operators. As 



Although fermions and hard-core bosons look alike 
in real space, in terms of their local density, one ex- 
pect their momentum distribution to be rather differ- 
ent from one another. Indeed the ground-state (GS) 
of a Fermi gas is the Fermi sea, in which all eigen- 
states - momentum states in a free gas - are filled 
up to the Fermi energy, whereas interacting bosons ex- 
hibit a peak in their the momentum distribution, around 
fc = 0. Indeed for hard-core bosons at zero-temperature, 
{■4)^x)^{x')) ^ \x — x'\^^l'^ when \x — x'\ — >• oo,'^^ and 
n(k) = / dxdx' e^^'^^~'^'^ {'^'^ {x)^(x' )) ^ /c^^/^ around 
fc = 0, as schematized in Fig. 3. This algebraic diver- 
gence is a result of the enhancement of quantum fluctua- 
tions in reduced dimensions, that forbid the existence of 
a true condensate. 



B. Free- fermions on a 2-leg ladder: exact solution 

We now turn to the two- leg ladder model Eq. (1.1) in 
the case of fermionic particles. Single-particle eigenstates 
of this Hamiltonian defined on a ladder with 2L sites are 
plane waves of the form: 



|fc,fc,) = 1 ^e-'=^-(|j,l)+e^^-''b-,2)) , (2.2) 
i=i 

where /c is a multiple of 27r/L, ky = 0, tt, and \j,£) is a 
state with a single particle at site The associated 

eigenvalues read: 

£fey(fc) = -2i cos(A;) - tj_ cos(fcy) , (2.3) 
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FIG. 4: (Color online) Left: single particle dispersion bands 
eo{k) (black curves) and e,r(fc) (red dashed) plotted (in units 
of t) for the two-leg ladder model Eq. (1.1) with t±/t = 1 (a) 
and t±/t — 3 (b). Right: Corresponding particle filling (per 
rung) {uj} — 2p plotted against the chemical potential n/t. 

which corresponds to the band dispersions depicted on 
Fig. 4. For fermions at zero-temperature, the average 
filling p vs chemical potential p is simply the sum of 
the fillings of each band. Two examples for t± = t and 
t± = 3i are given in Fig. 4. They correspond to two 
situations. First, if t± < 2t, the system has either four 
Fermi points (when — 2t + t± < jj. < 2t — t±) leading to 
a two-mode LL, or two Fermi points, leading to a single- 
mode LL equivalent to a chain. Cusps in the p{p) curve 
of Fig. 4 signal the change in the number of Fermi points. 

Secondly, when t± > 2t, the two bands are separated 
by a gap 

A = 2i_L - 4t , (2.4) 

and, at half-filling p = 1/2, the system is a band insula- 
tor. The corresponding GS has one particle on each rung 
in the symmetric state (|j, 1) -I- \j, 2))/\/2 with a total en- 
ergy Eqs = —Lt±. Away from the half- filling, the ladder 
is in a single-mode LL. 

The situation for hard-core bosons cannot be sim- 
ply understood in terms of filling single-particle states. 
A for a single chain, a Jordan- Wigner transformation 
can be performed by choosing a particular path for 
the string along the lattice. ''"^"''^ However, the fermionic 
model thus obtained cannot be solved exactly as it con- 
tains correlated-like hopping terms. Mean-field approx- 
imations yield reasonable but non-quantitative results, 
underlying the fact that the bosonic version of the model 
is non-trivial. In the following, we derive quantitatively 



exact results working in the bosonic language. 

III. ANALYTICAL RESULTS FOR HARD-CORE 
BOSONS 

A. Strong coupling limit 

1. Gap at half-filling 

In this section, we compute the gap at half-filling in 
the \aige-t± limit, in which t± corresponds to the en- 
ergy cost to add one extra hard-core boson on top of 
a half-filled ladder. Details of the calculation can be 
found in appendix A. Setting t = 0, rungs decouple 
and four states are available on each of them: an empty 
state \0)j with energy Eq = 0, two 1-particle states 
\l±)j = ^ (|j, I) ± \ j, 2)) with energies Ei± = Tt± ~ 

and a 2-particle state \2)j with energy E2 = —2p. These 
energies are plotted against the chemical potential /i in 
Fig. 5. The average filling p of the lattice is as long as 
p < — 1 for —t± < p <t± and 2 as soon as p > t±. 

We start from a half-filled ladder of hard-core bosons 
- p = and N = L bosons - and treat Hi = 
ei^l j^e j+i + ^^■'^■^ ^ perturbation of Ho — 

—t±J2j['^i,j^2.j + h-c.]. The GS of Hq is constructed 
by putting one particle on each rung in the symmet- 
ric state = ^ (|j, 1) + \j,2)). We have jiW) = 

|1+, 1+, . . . , 1+). Hi creates 2L particle-hole excitations 
on neighboring rungs and induces a second-order correc- 
tion to the GS energy: 

El = -tj_L~t^L/tj_. (3.1) 

The first excited states of Hq with + 1 particles are 
L times degenerate and best written in the momentum 
representation: 

u) 

(3.2) 

Hi lifts the degeneracy : indeed, it moves the extra par- 
ticle along the ladder with both nearest (first order per- 
turbation) and next nearest (second order perturbation) 
hopping and creates 2{L — 2) particle-hole excitations on 
nearest-neighboring rungs. Therefore the corrected ener- 
gies i?L+i(fc) of the L -\- 1-particle states are 

Ei+iik) = -tj_{L-l)-2tcos{k) 

- it^/tj_) cos(2A:) - (i - 2)tVi-L- (3-3) 

The situation and calculation are similar when doping 
the half-filled ladder with holes. The size of the plateau 
at half-filling is then 

A, = 2[^i+i(0) - El] = 2t^ -4t + 2t^/t^. (3.4) 

The gap is larger in the case of hard-core bosons than 
it is in the case of spinless fermions. To conclude this 
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FIG. 5: (Color online) Energies of the four available states on 
a given rung, in the limit t/tx — >■ 0. 



section, it is instructive to look at the perturbative ap- 
proach in the fermionic case, for which we know that the 
gap is exactly = 2t± — At, without a second order 
correction. This is due to the fact that the half-filled 
GS \L) = ■ ■ • L-i |0) and the L + f-particle states 
\{L + l)k) ~ ^ X]j are exact eigenstates 

of the full Hamiltonian. Hi does not create particle- 
holes excitations, since j = 0- the contrary, 
in the bosonic case, b'''^ jb^_^j = K^i.j^ij + + 
j^i j)- The diagonal terms vanishes because 
of the hard-core constraint, while the off-diagonal terms 
are equal because of the commutation of the bosonic 
operators (they vanish for fermions!). This difference, 
c\ ,c\_ ■ — versus b\ :b\ ■ = b\ -bl -, illustrates how 
the simplest quasi-lD system is indeed sensitive to quan- 
tum statistics. 



2. Effective model for incommensurate fillings 

We can analyze the strong rung coupling limit in the 
same spirit as for Heisenberg spin ladders in a magnetic 
field. ^■^'^^ According to Fig. 5, at a critical value of the 
chemical potential fic = t±, levels and \2)j be- 

come degenerate, with an energy Eq = —2t±. Doping 
the plateau with particles can be studied through degen- 
erate perturbation theory on the Hamiltonian 



Hi = -t^ Nj&£j+i+h-c.]-(Ai-/^c) 



'bl A 



(3.5) 



with respect to the GS Hamiltonian 

h.c] — IJ.C 



'^0 ^ -t^TAbl^b,^^ 



j 



3.1 



Let us call Pq the projector on the GS subspace, in which 
and \2)j are the only states allowed on a given 



rung j. We call Qq the projector on the complementary 
subspace. The effective Hamiltonian to second order in 
perturbation theory is'*'' 

Heff = PoHiPo + PoniQo[E - Hor^Qo-HiPo, (3.7) 

with E the eigenvalue of the degenerate subspace under 
consideration. Here, given the degenerate GS and the 
form of "Hi, a virtual excited state is a state with an 
empty rung. Then, in everything that follows, Qo[E — 
Ho]-'Qo = Qo[-2t^L + 2t^{L - l)]-^Qo = ~Qo/{2t^). 
Again, a more detailed calculation is given in appendix A. 
We find the following effective Hamiltonian expressed in 
terms of spinless fermions: 



"Heff = —t 



E B^..t 



h.c 



{P- - lJ'c)^{nj +1) 



3 



t' 

2f[ 



^2 



(3.8) 



In this language, a spinless fermion stands for a dou- 
bly occupied rung, while an empty site indicates a rung 
with one particle in the symmetric state. Note that 
one recovers the expression of the gap to second or- 
der. One particle is added when the effective chemi- 
cal potential equates the bottom of the energy band: 
^ — fj,c — 2t^/t±_ = —2t — t^/tj\_. On the other side 
of the plateau the exact same Hamiltonian is obtained, 
with spinless fermions now standing for empty rungs, and 
A^c = -t±. 

As already noticed in the context of 2D supersolid 
phases of quantum antifcrromagnets in a field'**^^*' or 
Heisenberg ladders in a field, '^^ an important ingredient 
which emerges as a second order process in the above 
effective Hamiltonian Eq. (3.8) is the correlated (or as- 
sisted) hopping between second neighbors. For Bose- 
Hubbard chains in the limit of large but finite on-site 
repulsion, an effective model of spinless fermions similar 
to Eq. (3.8) was obtained.^* We discuss further the rela- 
tionship with Bose-Hubbard chains and spin ladders in 
an external field below in sections V C 1 and VI C. 



B. Weakly coupled chains: bosonization 

In this section, we develop a low energy approach in or- 
der to analyze the behavior of two weakly coupled chains 
of interacting bosons at half filling. Our starting point is 
the following Hamiltonian % = Ho + H± with 

U 



(3.6) no = -tJ2 [bh+A,3 + ^■^■] + V E "^.j- ("^^^ - 1) 



3,e 



3-i 



Hi_ = -i±E [^I,.-^2,,+h.c.] 



(3.9) 
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T-lo describes two uncoupled chains of bosons interacting 
via the onsite repulsion U such that Eq. (1.1) is recovered 
for U = oo. couples the two chains. The low-energy 
excitations of a single chain of interacting bosons are col- 
lective excitations corresponding to the sound modes in 
the quasi-condensate. A suitable approach to describe 
them is bosonization.^ Within this representation, we can 
write a bosonic creation operator ?A^(a;) {ipl{x) being the 
continuum limit version of the lattice boson creation op- 
erator b^ ) as 



-.1/2 



p v$ 

TT 



gi2p(irpa;-*(a;))g-ie(a:)^ 



(3.10) 

where p is an integer and p is the boson density. $ and 
9 are two dual bosonic fields obeying the following com- 
mutation relation 



$(a;), -V6'(a;') ^i5{x-x'). 

TT J 



(3.11) 



In this language, — - V$ accounts for the long wavelength 
density fluctuations, while 9 can be regarded as the su- 
perfluid phase. Using Eq. (3.10) and taking the contin- 
uum limit, we can rewrite Hq in Eq. (3.9) as^ 



Ho 



V 

2^ 



E 



dx 



K{V9i 



K 



(V$, 



(3.12) 



where we have introduced the pair of bosonic fields 
and 9^ living on chain I. The low-energy physics of each 
chain is characterized by the same two parameters, v the 
sound velocity and K the Luttinger parameter which en- 
codes interactions in the system. For interacting bosons 
with on-site repulsion, K > \. With nearest-neighbor 
interactions it is possible to reach K thus making a 
connection with XXZ spin chains. Using Eq. (3.10) we 
can rewrite Hx^ as 



Hi 



dx cos( 



(3.13) 



1 + 4 cos(27rpa; - <i>2 - $i) cos($2 - <I>i) 

-2 cos(47rpa; - 2$2 - 2$i) + 2 cos(2$2 - 2$i) 



In order to write this equation, we restricted the sum- 
mation in Eq. (3.10) to the values p = 0,1,-1 which 
provide the leading terms in the perturbation. Half- 
fiUing corresponds to p = l/(2a) with a the lattice spac- 
ing. Introducing symmetric and antisymmetric fields 
defined as = ($i + ^2)/V^:^a = (^'i - $2)/\/2, 
9s = {9i + 6'2)/v/2, 9a = (6*1 - 92)/\/2 and retaining only 
the non-oscillatory terms, we can rewrite T-L = "Ho + "Hj^ 



E 

ttq; 



^ / dx 



dx 



cos 2 COs(\/29a) COS (%/8$^) 

fzcos (V20,)cos(V8$,) , (3.14) 



where Ka = Kg = K and Va = Vg = v. We see that 
the coupling between the two chains is quite complex 
in the bosonized representation. To study this kind of 
Hamiltonian a strategy consists in deriving RG equations 
for the various coupling constants. To do so, we first 
write the Euclidian action S associated to Eq. (3.14) in 
the following general form: 



S 



E 

i=s,a 



1 



27rif, 



dxdr 



1 



dxd{vT) 
2ya 

+^ 

Tra^ 



{V29a) 

/I U," 

COS {V29a) C0s(V8$a) 

cos 

■ cos(\/8<I)s 



(3.15) 



where we have introduced four dimensionless coupling 
constants y±, yg, ya and finally which is generated 
under the RG flow as we will see. These couplings con- 
stants have bare values which are directly extracted from 
the Hamiltonian in Eq. (3.14) and read yA_{a) = at±/v, 
2/s(a) — ya{a) — at±v, gsia) = 0. We write RG equations 
for the six coupling constants using standard techniques^: 



dy±_ 

dl 
dya 
dl 
dys 
dl 
dgs 

dl 

dKa 
dl 

dKg 
dl 



1 



2Ka 

1 
1 

2K'a 



y±ii), 

-2Ka)ya{l). 



(3.16) 



= (2~2Ks)gS) + -y^{l)yg{l), 

yl{i) , ym 

4^4' 
= -AKl{l)yl{l)~2gl{l)Kl{l), 



The running short distance cutoff is parametrized as 



a{l) 



Since the bare values of Kg a are larger than 



one, the coupling y±^ is strongly relevant compared to 
the other couplings and y±^ is therefore driven to strong 
coupling. This implies that the fleld 9a becomes gapped. 
The superfluid phase fields of both chains lock together 
as soon as interchain tunneling is switched on. Next we 
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introduce the scale l± such that y±{l±) = 1 which equiva- 
lently defines an energy scale ~ . For I < l±, one 
can use the full set of RG equations written in Eq. (3.16). 
For I > l±, one can simplify the action (3.15) by replac- 
ing the field cos{^/29a) by its average value (a similar 
procedure is described in Ref.52 for coupled spin chains) 
and thus obtain an effective action valid at lower energy 
scales I > l±. The phase field 9a being gapped, we focus 
on the symmetric (s) sector whose effective action 
now takes the simpler form 



eff 



dxdTi - 



1 



'9s_ 

7ra^ 



1 

Vs 
COS 



(a,$,)2 + «,(a,$,)2 

iVs^s)} , (3.17) 



where we have introduced Kg = Ks{li_) and gs — gs{l±)+ 
'2r]ys{l±) with r; = {cos{V29a)) = 0(1). The action Sf 
is of the sine-Gordon type. The RG equations associated 
to Kg and Ijs are obtained by further integrating high- 
energy degrees of freedom between l± and I, leading to 



dl 

dK, 
dl 



(3.18) 
(3.19) 



which are the standard Kosterlitz-Thouless flow equa- 
tions.'' The coupling gs flows to strong coupling when 
2Ks{l) is smaller than 1. Since Ks{l) always decreases 
during the fist stage of the RG transformation, gs is 
driven to strong coupling and, ultimately, a gap opens 
up in the symmetric sector, at a very low energy scale. 
We stress that the gap opening is non trivial from the 
RG point of view: first, a gap opens in the antisym- 
metric sector (pining of the field 9a)-, which eventually 
induces the opening of the gap in the symmetric sector 
(pining of the field <i>s). The scaling of the gap with t±, 
obtained from the numerical solution of the RG flow, is 
shown in Fig. 6. It strongly depends on the initial value 
of K. For hard-core bosons ~ K = 1 - the charge gap 
grows exponentially slowly with t±: 



(X e 



(3.20) 



Indeed, during the first step of the flow, I < l±, ys is 
irrelevant and renormalized downwards to a small value, 
while Kg decreases to a value smaller but close to 1. At 
the start of the second step, ys has become a relevant 
perturbation but is still very close to marginality. On 
the contrary we could imagine a situation with K well 
below 1 - for instance, adding nearest neighbor repulsive 
interactions in (3.9) - and j/s would much farther away 
from marginality at the start of the second step, therefore 
leading to a much faster - power-law - opening of the gap 
(see Fig. 6), 



K = I 
exponential fit (3.20) 




□ K = 3/4 




power-law fit (3.21) 


P ■ 
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□ e' 
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FIG. 6: (Color online) Charge gap obtained from the numer- 
ical integration of the two-step RG flow. The scaling of the 
gap depends strongly on the initial value of K. We find for 
the exponent of the power-law, eq. (3.21), h = 2.81 in good 
agreement with the DMRG result of equation (4.9). 



Finally, this whole RG analysis shows that a gap in the 
symmetric sector opens up as soon as transverse hopping 
is switched on. This result is drastically different from 
the fermionic case where a band gap appears only for a 
finite large value of 



C. Comparison with spin- 1/2 ladders 

A two-leg ladder of hard-core bosons is equivalent 
to a spin-1/2 two-leg ladder. Indeed, spins-1/2 satisfy 
mixed commutation relations, when expressed in terms 



of S. 



^ ox _|_ ■ qV 



A, CX {t^/tf 



(3.21) 



A hard-core boson on a given 
site corresponds to a magnetic moment with = -1-1/2 
while an empty site corresponds to the opposite with 
= —1/2. Consider the model of hard-core bosons we 
have studied and add nearest-neighbour repulsion along 
and between the chain. It maps to a spin-1/2 ladder: 

L 

For = = it reduces to the model Eq. (1.1), 
with t — J/2 and t± = J±/2. The strong coupling 
analysis^'^^^^'^"^ leads to a similar picture: singlets form 
on each rung and a finite spin gap (m = magnetization 
plateau) extends around zero magnetic field. A weak- 
coupling analysis leads to the same conclusion. However 
the mechanism for the formation of the plateau is dif- 
ferent than in our study of hard-core bosons because of 
the exchange coupling. We recall below the bosonized 
form of the Hamiltonian (3.22). In a similar way to the 
physics described in section IIIB, the dynamics decouple 
into two modes, symmetric and antisymmetric such that 
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one can write H = 'Hs + Ha with 



2tt 



dx 



(2iray J 



(3.23) 



'Ha — 



27r 



dx 



KaiVea)^ + --iV<i>af 



2ga 

(27ra)2 

"(27ra)2 



J dx cos[\/8$a] 



da; cos[-\/80a] 



(3.24) 



Here, = K [l ^ ^) , = K (l + ^) , = 

f/a = aJ^,Ga = T:J±a. K is the original Luttingcr pa- 
rameter of the single spin chain and depends on , the 
anisotropy parameter. For instance K — \ for = 
(hard-core bosons) and iiT = 1/2 for = J (Heisenberg 
chain). Note that the z exchange along the chains should 
also give rise to an umklapp term of the form: 

Uu = / rfxcos[\/8$,]cos[y8$a] (3.25) 

(27ra)2 J 

with gu (X J'^ . However as long as K > 1/2 this term is 
always less relevant than the other interaction terms and 
will not open any gap. When vanishes one should also 
include the term cos[\/29a] cos[2a/2$s] which we studied 
in the previous section. For a non-zero J^, it is anyway 
less relevant and the gap in the symmetric sector opens 
directly because of the exchange along z. Indeed, for all 
< 1 and ^ 0, Kg < 1 and gs is relevant, therefore 
ordering the field $s. The scaling of the gap is known 
from the sine-Gordon model and depends on the initial 
conditions for Ks and gs- Away from marginality, the 
charge gap scales as a power-law: 



A, 



(3.26) 



Carlo^'' (QMC) and density-matrix renormalization 
group6""^2 (DMRG) results at half-filling. The Rung- 
Mott insulating state (RMI) will be studied in great de- 
tails using these two techniques. 



A. Density plateau at half-filling 

1. QMC results 

We use the QMC stochastic series expansion (SSE) 
algorithm,'"'^ in its directed loop framework. ^^'^'^ This 
method has been used quite intensively for the last 
decade. For our purpose here, the implementation of 
the algorithm is rather straightforward since one can, for 
simplicity, work in the representation of the equivalent 
spin-1/2 XY model in a transverse field. Indeed, using 
the Matsubara-Matsuda mapping, ^"^ hard-core bosons 
operators are replaced by spin-^ operators: 



-2t^ {Sl^Sl^ + Sf jSl^) + const. (4.1) 



This equivalence with spin-1/2 systems, already evoked 
in section HI C to discuss the connection with XXZ spin 
ladders, will also be used later in section VI when dis- 
cussing potential links of our results with experiments 
on magnetic field induced LL phases in spin ladders. 

Coming back to the bosonic language, we work in the 
grand-canonical ensemble where the total particle num- 
ber (longitudinal magnetization) can fluctuate. We mea- 
sure the following observables: (i) The density of particles 



2L 



(4.2) 



2L 



3/ 



For instance, at the isotropic Heisenberg point J^^ = J, ^--j-j compressibility 
K = 1/2 and the gap opens linearly with 

Since the antisymmetric sector is always gapped, we 
expect the familiar picture of "rung-singlet" known for 
Heisenberg ladders^^ to hold here for the generic U(l) 
symmetric case for which the gapped state is called 
a Rung -Mott- Insulator (RMI). However, the correlation 
length of the symmetric mode ^ ~ l/A^ is much larger 
than for Heisenberg ladders. Therefore, the usual short- 
range resonating singlets picture, very useful for spin lad- 
j^gj.g 55,56 jjg^g |-p taken more carefully here, especially 
in the limit of small t_\_/t. 



J2 



(4.3) 



and (iii) the superfluid stiffness 



1 d'Eoi^) 
2L (9^2 



(4.4) 



ip=0 



IV. RUNG-MOTT INSULATOR: NUMERICAL 
RESULTS 

Coming back to the original problem of hard-core 
bosons, we present in this section quantum Monte 



In the above definitions, /3 is simply the inverse tem- 
perature fceT, and (/? is a small twist angle enforced on 
all bonds in the direction^^'®^ along the legs. Techni- 
cally, the superfluid stiffness^^ is efficiently measured via 
the fluctuations of the winding number^^ during the SSE 
simulation. In Fig. 7, we show QMC results for the to- 
tal particle density p, the superfluid density psi (directly 
related to Tgf as discussed below in Sec. IV B), and the 
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FIG. 7: (Color online) QMC results for the particle density 
p (top), the SF density psi (middle), and the compressibil- 
ity K (bottom) obtained on 2 x 64 ladders with t± = t = 1 
(Left) and t± = 2t (Right), for two inverse temperature 
j3 = 128, 256, as indicated on the plot. On the left panel 
{t± = t), the compressibility data are obtained from either a 
direct measurement of k [full lines Eq. (4.3)], or from a nu- 
merical derivative of p{fi) [circles Eq. (4.5)]. Note for t_L = t 
(left) the superfluid density data are also shown for L = 128 
and L = 256 at /? = 4 X L. 



compressibility k. Note that the later can also be ex- 
tracted from the numerical derivative 



k(^) 



dp 
dfi' 



(4.5) 



QMC data are shown for two representative values of the 
interchain hopping: t± = t and t± = 2t. As expected for 
large t±, a plateau in the density at half- filling is clearly 
visible for t± = 2t if —0.5 ^ p./t ^ 0.5, signaling the 
incompressible RMI state where the superfluid density 
and the compressibility both vanish. For t± — t, inter- 
pretation of finite size numerical data is more difficult 
since the gap at half-filling turns out to be much smaller. 
Indeed, the density curve versus ^ does not show any 
visible plateau (top left of Fig. 7). Nevertheless, the gap 
tends to show up in the compressibility k which displays 
a downward feature, signaling the insulating state. The 
SF density also vanishes, as it becomes visible when L is 
increased from 64 up to 256. Obviously, when t±/t gets 
smaller, it will be quite difficult to directly get quantita- 
tive information from p or k, like for instance the value 
of the gap As at half-filling. If we want to do so, one 
needs to run QMC simulations in the GS, i.e. at temper- 
ature well below the energy gap T <^ and in the ther- 
modynamic limit, i.e. for system lengths L 3> ^ where 
^ ~ 1/As is the correlation length associated with the 
short-range order of the RMI. Of course this would be a 
difficult task, that we can fortunately circumvent. 

In a QMC framework, the most useful observable to 
characterize the insulating state, and to get a precise es- 
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FIG. 8: (Color online) Charge gap of the two-leg ladder 
model as a function of the interchain coupling obtained from 
DMRG calculations, and compared with perturbation theory 
Eq. (2.4) for free-fermions and Eq. (3.4) for the strong rung 
coupling expansion. The red curve is a fit to the expression 
Eq (4,8) (see the text). 



timate of the correlation length ^ is the superfluid stifl- 
ness computed directly at half-filling, which is expected 
to vanish as exp(— i/^). A direct estimate of the gap with 
SSE, in principle possible, would require a much more im- 
portant numerical effort. Zero-temperature DMRG sim- 
ulations are more competitive to measure tiny gaps as we 
discuss now before coming back to SSE estimates of the 
SF stiffness and the correlation length. 



2. Gap from DMRG 

The gap is computed at zero temperature with DMRG 
using the definition 

A, ^ EoiN + I) + Ea{N - I) - 2EoiN) , (4.6) 

where Eq{N) is the GS energy with N bosons, so that A^ 
is directly the width of the plateau. In the DMRG cal- 
culation, we keep 1000 states per block and extrapolate 
the results over sizes ranging from 24 to 160 using the 
following ansatz from the finite-size scaling of the gap 



As{L) = A, 



A 

r 



(4.7) 



The result is reported on Fig. 8 in which the strong- 
coupling of Sec. Ill A 1 and free-fermions predictions are 
also given for comparison. According to the bosonization 
prediction of Sec. IIIB, we fit the opening for t±^ < 2t 
using the following law 



A. = A" 



exp 



(4.8) 
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FIG. 9: (Color online) Comparison of the opening of the 
gap in the XXZ ladder for three different situations. DMRG 
results (symbols) are compared to Eq. (4.8) for hard-core 
bosons, to Eq. (4.9) for coupled XXZ chains, and to the linear 
behavior from Greven et alJ''^ in the SU(2) case. 




FIG. 10; (Color online) SF density p^f at half-filling (normal- 
ized by the decoupled chains value plf — I/tt), plotted versus 
t± . QMC results obtained for ladders of size 2 x L in the GS 
at inverse temperatures /3 = 41/ with L = 32, ■ • ■ , 512. 



B. Superfluid density 



with Ag, a and as fitting parameters. The obtained 
values are respectively ~ 15.624i, a ~ 5.294 and f]^ ~ 
0.0075i. The latter critical value is perfectly compatible 
with = within our numerical precision. The opening 
of the gap is particularly slow (non-analytic in t±), but 
reaches a sizable magnitude 0.074t ^ A^ ^ t for t < t_L < 
2t, thus providing a first clear quantitative difference with 
free fermions. 

Since it is well known that SU(2) spin-1/2 ladders have 
a gap at the isotropic poinf^'' J = J_l, and that this gap 
opens linearly^" with J± in the limit of weakly coupled 
chains, we investigate how interactions modify the open- 
ing of the gap in order to have an intermediate situa- 
tion. We introduce intra-chain interactions with nearest- 
neighbor density interactions corresponding to = 0.5 
in the XXZ language model. From the Bethe-ansatz solu- 
tion, the Luttinger parameter of the chains at zero inter- 
chain coupling is thus K = 3/4. We also recompute the 
SU(2) (with .F ^ Jl = l and for which K ^1/2) gap to 
check the numerical accuracy and gather results in Fig. 9. 
As expected from the RG calculation of Sec. Ill B, the in- 
termediate opening is best fitted by a power-law with an 
exponent of the order of three: 



0.2234 t 



(4.9) 



close to the exponent obtained from the RG calculation, 
and non-trivially related to the value of K (see Fig. 6). 
The gap in the isotropic limit reproduces well the earlier 
findings of Ref.70. The law governing the opening of 
the gap is thus very sensitive to interactions, with their 
expected tendency to boost the charge gap. 



1. QMC results 

As introduced by Fisher, Barber and Jasnow in Ref.67, 
the superfluid stiffness (or helicity modulus) Tsf is de- 
fined by imposing twisted boundary conditions in one 
direction, such that the hard-core boson Hamiltonian 
Eq. (1.1) reads now 



(4.10) 



Such a twist angle tp mimics the effect of a phase gradi- 
ent imposed in the bosonic wave function. In a superfluid 
state, it would lead to a superflow with a gain in the ki- 
netic energy density'^'^ SEo{f) = jT^f-^^ = ^psf(p^, 
where m* is the effective mass of the bosons, m* is de- 
termined in the diluted limit where the cosine dispersion 
2t cos fcj. can be approximatively described as a quadratic 



h-'kt 



free bosons dispersion , 
therefore 



thus leading to 



Psf = 



1 1 d^Epjip) 
2L 2t difi^ 



ip=0 



i, and 



(4.11) 



Based on the winding number fluctuations, ' we com- 
pute for various system sizes the superfluid density at 
half- filling fi — versus interchain hoppings t±. As al- 
ready discussed, it is crucial to get zero-temperature es- 
timates. This is ensured here by performing SSE simula- 
tions at /3 = 4 X L for the values of t± considered in the 
following analysis. In Fig. 10, the SF density psf is plot- 
ted against tj_ for L = 32, 48, 64, 96, 128, 192, 256, 512, 
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FIG. 11: (Color online) (a) SF density plotted versus the 
linear size L for various values of t^, as indicated below with 
different symbols. The linear-log plot clearly shows an expo- 
nential decay with L. (b) Collapse of the above data obtained 
by rescaling the x-axis L ^ L/^ in order to get all points on 
a single curve. The full line is cx exp(— L/^). The obtained 
correlation length ^{t±) are displayed in Fig. 12. 
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FIG. 12: (Color online) Correlation length of the bosonic 
ladder plotted versus the inter-chain hopping t_i_. The points 
obtained from the data collapse of Fig. 11 (black circles) are 
well-described by the exponential form displayed on the plot 
(green curve). 



scaling regimes are clearly visible. With this technique, 
we have been able to extract the value of ^ as a function 
of t± down to t± = OAt. For smaller values, it was almost 
impossible to see any downward curvature. 



and normalized by its tj^ = value which is simply, in 
the thermodynamic limit^^ plf — 1/tt. We observe in 
Fig. 10 a slow decay of the SF density when the inter- 
chain hopping is increased. Even for the largest system 
considered here (1024 sites), the SF response is still very 
important for t± < t/2. While for large enough values 
of t± ~ t, the gapped nature of the system is obvious, 
with a vanishing SF density, it appears quite difficult to 
predict the existence of any critical coupling f]^ from the 
actual behavior of Psi{t±), as displayed in Fig. 10. 



2. Scaling analysis 

We therefore use another strategy to extract the criti- 
cal properties. As visible in Fig. 11(a) for various values 
of t±, the SF density decay exponentially at large L, at 
least for not too small values of the coupling: t± > 0.7t, 
while it roughly remains constant for the sizes considered 
here when t± < 0.7t. We therefore identify two scal- 
ing regimes, depending on the scaling variable x = L/^, 
where ^ is the correlation length of the RMl state: 

^-/^-' - { exp(-x) (4.12) 

Data of Fig. 11(a) have been analyzed by rescaling the 
X-axis L — )■ L/^ in order to collapse the data : using the 
best fit of the form Aexp{—L/^), we have fixed ^ = 61 
for t± = t and adjusted all the other data to get the best 
collapse. Results are displayed in Fig 11(b) where the two 



C. Correlation length 

In order to convince oneself that the data collapse 
technique gives rather good estimates of the correla- 
tion length, we give in Appendix D a benchmark of the 
method in the case of the integrable t — V model. Now 
we show how the correlation length of the bosonic lad- 
der, extracted in the data collapse of Fig. 11, varies with 
t±. We report the estimates in table 1, and plot In^ 
against t± in Fig. 12 where two fitting functions of the 
form ^oexp [a,t/ {tj_ — t']_)] are displayed. When leaving 
free the three parameters ^o, a, and t]^, we obtain for 
the best fit (green curve in Fig. 12) a very small critical 
hopping ^5. = 0.007t, obviously compatible with the KG 
result f]^ — 0. When forcing t"]^ = in the fitting pro- 
cess, we also obtain a very good agreement (dashed blue 
curve) , which gives a very good support to the RG result 

e = 6exp(^a^). (4.13) 

Note that we can compare ^ with the inverse gap l/A^ 
extracted from DMRG simulations (Eq. (4.8) and Fig. 8). 
They both agree very well, showing in particular a rather 
large value of the parameter a ~ 5.2(1). As already em- 
phasized, the physics of short-range correlations - the 
usual hallmark of quantum spin ladders - is not real- 
ized here, except in the limit of large interchain hopping 
tj_ ^ t. It is indeed quite instructive to compare the 
results of ^ for the bosonic ladder model with the cor- 
relation length of SU(2) symmetric Heisenberg ladders 
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t± 
t 


[bosonic ladder] 


[SU(2) ladder] 


1 


60(1) 


1.99 


0.9 


106(2) 


2.28 


0.8 


2.15(5) X 10^ 


2.65 


0.7 


5.9(1) X 10^ 


3.12 


0.6 


1.9(1) X 10^ 


3.76 


0.5 


1.2(1) X 10* 


4.64 


0.4 


1.6(1) X 10^ 


5.93 



TABLE I: QMC estimates for the correlation length ^ of the 
half-filled RMI. These data are also plotted against in 
Fig. 12. For comparison, the third column shows the correla- 
tion length of the corresponding rung singlet state of SU(2) 
Heisenberg ladders [model Eq. (3.22)] with J± = = t± 
and = J = t (see also Fig. 9). 



(table I), known to display very short-range correlations, 
thus providing the useful picture of short-range resonat- 
ing valence bonds. In the bosonic ladder case, the RMI 
state has correlations along the legs which extend over a 
very large scale ^, even for t± ^ t, as one can read from 
table I. Practically speaking, for instance in an optical 
lattice setup made of at most several hundreds sites, two 
weakly coupled chains will behave as a gapless system 
(see the discussion in section VI C) . 



V. LUTTINGER LIQUID BEHAVIOR 

We now turn to the case of incommensurate fillings, 
and study the behavior of the Luttinger parameter Kg 
of the symmetric gapless mode which governs all corre- 
lators. It has a non-trivial dependence on density and 
interchain hopping which requires non-perturbative nu- 
merical estimates : two complementary approaches are 
used from DMRG and QMC calculations. Analytical re- 
sults for several limiting cases are also presented. A first 
account of this behavior was given in Ref.24. We improve 
the numerical accuracy of the calculation and we present 
a detailed discussion of the physical origin of the curves. 



Luttinger parameter from the correlation 
functions: DMRG results 



In this section, we describe how Ks is extracted from 
the computation of correlation functions using DMRG 
with open boundary conditions (OBC). As Kg governs 
the decay of all algebraic correlators, the idea is to find 
a suitable order parameter which provides a direct ac- 
cess to Kg. The most natural one would be the bosonic 
propagator along a chain {b\ j^^b^ j), which has a lead- 
ing algebraic decay at zero-momentum governed by the 
exponent 4:/ Kg. However, numerical calculations shows 
that the signal is spoiled by subleading oscillating terms. 
A better choice is actually to look at the propagator of a 
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FIG. 13: (Color online) Typical rung-pairing correlation func- 
tion P{x) computed with DMRG and fitted by Eq. (5.2). The 
effect of the interchain hopping t± is strong on both the am- 
plitude and the exponent. 



pair of bosons created on a rung, defining: 



(5.1) 



This correlator has a zero-momentum leading term with 
a decay exponent I /Kg and do not display sizable sub- 
leading corrections as shown in Fig. 13. The data are 
here calculated from the middle to the edge of a system 
with L = 192 and keeping 2000 states. With OBC, the 
wave-function naturally vanishes at the edges, causing a 
fall of the correlator when x,y ^ L. In order to get the 
best estimates for Ks, we have to take into account this 
finite-size corrections. Thanks to the conformal field the- 
ory results of Refs.24,72 and which already proved to give 
very good fits of DMRG data,^^'^'^'^''' we use the following 
fitting function for a superfluid-like correlator: 



Pix,y)^A 



y/d{2x\2L)d{2y\2L) 
d{x + y\2L)dix ~ y\2L) 



1 i/Ks 



(5.2) 



in which we take y = L/2. Here, d{x\L) — ^ sin (^) is 
the cord function^^^. Typical fits are displayed in Fig. 13, 
from which a rather accurate decay exponent can be ex- 
tracted. 

Using this procedure, we mapped out the behavior 
of Ks as a function of the density p for various inter- 
chain hopping t±. The behavior is plotted in Fig. 14 (in 
a symmetrical way because of particle-hole symmetry) 
and shows several striking features. Kg is strongly af- 
fected, as expected, by the vicinity of the commensurate- 
incommensurate (C-IC) transition occurring at p — 1/2. 
Between and 0.5, a strong asymmetry is found, which 
is surprisingly qualitatively very similar to the behavior 
of the Luttinger parameter of the charge mode in the 
ID fermionic Hubbard model, and for which we will 
give supporting arguments. The limiting value of Ks is 
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FIG. 14: (Color online) DMRG results for the Luttinger pa- 
rameter Ks plotted versus the density p. Various symbols 
represent different interchain hopping strengths t±/t. The full 
black line is the perturbative results Eq. (5.11) for t^/t — 8, 
as derived in Sec. V C 2. 

naturally one in the limit of uncoupled chains^^^. Be- 
fore embarking in a more careful analysis, we may notice 
that the C-IC scenario tells that the expected value of 
Ks approaching the commensurate lines by varying the 
density is half that of the critical Ks obtained by varying 
the coupling which does yield Ks = 1/2 in the vicinity of 
p=l/2. 

B. QMC estimate of the Luttinger parameter 

The Luttinger parameter of the gapless (symmetric) 
mode can also be computed quite efficiently using QMC. 
While the Green's functions are in principle accessible 
with SSEJ^ a precise determination of long distance off- 
diagonal correlators of the form '^'-'^ ^^^Y 
task. An interesting alternative using QMC is to take ad- 
vantage of the "hydrodynamic relation" of the LL where 
phase and density fluctuations are related through 

Ks = Try/ K.Ts{. (5.3) 

As discussed above, it is straightforward to compute 
the superfluid stiffness Tgf and the compressibility k. We 
therefore compute the effective LL parameter Kcs{T) — 
TT k{T)T sf{T) . Finite temperature results are shown 
for a 2 X 64 ladder with t± — t in Fig. 15 where one 
sees three different regimes, (i) Kcff{T) = at high tem- 
perature because the SF stiffness is zero, (ii) When T 
decreases Kce{T) increases and saturates as soon as T 
becomes of the order of the finite size gap A(L). (fii) 
For T below A(L), oscillations in the compressibility k 
start to appear, thus signaling the appearance of small 
magnetization steps in the GS of the finite length sys- 
tem. This effect naturally induces oscillations of Kcs{T). 




FIG. 15: (Color online) Effective LL parameter Kcb{T) = 
tt^k{T)Tsi{T) computed with SSE simulations on a 2 x 64 
ladder with t± — t = 1 and for various temperatures, as in- 
dicated on the plot. The compressibility is computed either 
through the numerical derivative dp/dfj, (a) or directly using 
the SSE estimate Eq. (4.3) (b). Note that large oscillations 
visible at the lowest temperatures disappear when approach- 
ing the transition points pc = ±3. 



Note that these oscillations are more pronounced when 
K is directly computed through Eq. (4.3) whereas tak- 
ing the numerical derivative dp/ dp, has a smoothing ef- 
fect, as visible in Fig. 15 (a-b). It is however interesting 
to note that close to the diluted limit, where the level 
spacing at the bottom of the cosine band becomes much 
smaller A(L) ~ with z = 2, these oscillating fea- 
tures tend to disappear, even at the lowest temperature 
considered here. It appears very clearly here that the 
Luttinger parameter tends to 0.5 in this diluted limit, as 
well as in the other limit, when the RMI is approached. 
Below, we give analytical support to such a limit. In 
between, the incommensurate regime is characterized by 
a /i— dependent Luttinger exponent Ks{p) which takes 
non-universal values in the range (0.5— 1). These results 
are in very good agreement with DMRG estimates dis- 
cussed above, as shown in Fig. 16 where QMC together 
with DMRG results are displayed for tj\_/t = 2. Again 
one sees that a very large inverse temperature /3 is nec- 
essary to capture the GS properties close to the diluted 
limit. In particular this is quite visible in Fig. 16 for the 
SF density p^f which displays important finite tempera- 
ture effects in this diluted regime. 
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FIG. 16: (Color online) QMC results obtained for tj_ = 2t 
and 1/ = 64 at various temperatures. The full lines show the 
total particle density p whereas the SF part p^f is shown by 
small circles for the three values of (3 using the same color as 
for p. The Luttinger parameter Ks, extracted using Eq. (5.3), 
is plotted also for the three considered temperatures (colored 
diamonds), and compared with DMRG estimate (magenta 
squares) . 



We now give some physical interpretation which sup- 
ports the numerically observed Luttinger parameter Kg 
in various limiting cases. 



C. Limiting cases: analytical results 

1. Dilute Bose gas limit 

In the extremely diluted limit, the bosonic density 
grows with the chemical potential ii as the density of 
free fermions would do upon filling up the bottom of the 
cosine 0-band Eq. (2.3). This is apparent in the srong 
transverse hopping limit, where the Hamiltonian is actu- 
ally mapped onto a band of spinless fermions in equation 
(3.8). At low densities, interactions and correlated hop- 
ping can be neglected, and fermions are indeed free. In 
the weak coupling limit, a look at bosonic propagators 
helps to draw a similar picture. Considering symmet- 
ric bj,+ = {bj,i + bj^2)/V^, and antisymmetric bosons 
— {bj i — bj_2)/V2, we have: 



cos 



(5.4) 



(^U^)^(0)>-(f)'^(cos(^ 
and 

(5t.(.)6_(0))^(^)^(sin(M^ 

^ (5.5) 

Since interchain hopping is a relevant perturbation and 
the field da is gapped, the average of cosines decay ex- 
ponetially to 1, over a length scale ^a, while the average of 



Oa{0) 



sines decay exponentially to over the same length scale. 
In the next paragraph, we discuss in detail the behavior 
of the important length scale with t± and the filling 
p. Here it is enough to realize that at very low fillings, 
only one gapless mode remains, consisting of a Luttinger 
liquid of symmetric bosons. In Fig. 17, we computed 
both the filling of hard-core bosons and free fermions on 
a two- leg ladder with t = t±. It appears that they be- 
come identical in the extremely diluted limit, confirming 
both our strong and weak coupling analysis. This leads 
to the following square-root behavior for the density 



pip) 



1 

2^ 



which yields a compressibility 



(5.6) 



(5.7) 



Also, for very diluted systems, it is natural to expect 
that the SF fraction will tend to unity^^ : psi/p ^ 1- 
Again this is exemplified in Fig. 17 where Psf/phcb 1 
as Phcb ^ 0- Therefore, in such a limit the stiffness 
Tgf — 2tpsf is simply given by 



t p~ He 



(5.8) 



which, using the hydrodynamic relation for the Luttinger 
parameter Eq. (5.3) simply becomes 



(5.9) 



This is a striking difference between free fermions and 
hard-core bosons on a two-leg ladder: the former have 
a single-mode LL with K — 1 even in the vicinity of 
p — 0,1. Note that symmetric bosons indeed behave as 
ID hard-core bosons at very low fillings since their prop- 
agator now reads {b\_{x)b^{0)) ~ From Fig. 17 
we also notice that as the density increases, the follow- 
ing inequality phcb > Psi > Ps is verified, implying that 
Ks > 1/2 as the system is filled up. 



2. Doped RMI: strong rung coupling 

Starting from the effective spinless fermions Hamilto- 
nian derived in Sec. HI A 2 to second order in t/tj_, equa- 
tion (3.8), we deal with a ID chain for which it is possi- 
ble to compute both the Luttinger parameter K and the 
renormalized velocity v using the compressibility and the 
phase stiffness. Indeed one can easily compute the GS en- 
ergy of the effective Hamiltonian, to first order in t/t± 
and its response to a change in the number of particles 
or a twist in the boundary conditions (see for instance 
Ref. 78 as well as appendix C). We find: 



K = 1 



2t sin(7r/5) 



(5.10) 
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FIG. 17: (Color online) Various densities plotted vs /^/t for 
t± — t: hard-core bosons (QMC, red circles), free-fermions 
(exact, dashed line), and the SF density for the hard-core 
bosons (QMC, open circles). QMC results obtained for 2 x 64 
ladders at P/t — 512. Inset: SF fraction psf/Phcb (open blue 
circles) and fermion-boson ratio pa/phch (red circles) shown 
close to the critical point ~ —3t. The lines are guides to 
the eyes. 



with p = 2{p — 0.5) the density of fermions, that is, the 
density of doubly oecupied rungs. We now turn to the 
question of the relation between K and in the origi- 
nal system. As argued earlier, intuition suggests that the 
fermions in the effective model are related to the symmet- 
ric bosons of the original ladder, b+ = {bi +b2)/V^, and, 
from the related propagator (5.4), we saw that at large 
enough distances, x ^a, these particles actually behave 
as a Luttinger liquid. In the strong coupling regime, us- 
ing the mappings of Appendix B, we calculate the prop- 
agator to be {b^^{x)b^{0)) ~ x~^^^'^^\ We can readily 
make the identification K — 2Ks in the strong coupling 
regime and obtain the following perturbative expression: 



1 



t sin[27r(p - 0.5)] 



(5.11) 



The latter expression seems already a good approxima- 
tion for t_L — 8t, as visible in Fig. 14. Note that a 
similar calculation was first carried out by Cazalilla^^ 
in a different context. There, the author computed the 
Luttinger parameter of the Bose-Hubbard model in the 
limit of strong on-site interaction by mapping the prob- 
lem onto a chain of spinless fermions with attractive 
nearest-neighbor interactions and next-nearest neighbor 
correlated hopping. A parallel can be drawn in our case. 
One could argue that the symmetric bosons are effec- 
tively soft-core bosons, since one can put two of them on 
the same site. However it is important to note that one 
cannot put more than two, and that a perturbative cal- 
culation on the Bose-Hubbard model for fillings greater 
than one would lead to a different effective Hamiltonian. 

Naturally, by doping the rung-Mott insulator with 
holes, one recovers a symmetric expression for Ks when 




FIG. 18: (Color online) Ratio (,a/d as a function of the density 
p for different values of t±/t. 



p < 0.5: 



t sin[27r(0.5 - p)] 



(5.12) 



This being said, we now have a quantitative description 
of the behavior of Kg for large t±/t. 



3. Weak coupling limit and asymmetry of Ks 

We now provide a simple argument taken from the RG 
analysis to explain the qualitative behavior of Kg at low 
tj_/t, in particular the large asymmetry about quarter 
fillings. For t± < t, the length we introduced earlier 
can be computed from the RG: 



'2K/{4K-1) 



(5.13) 



where v — 2tsin(7rp) is the sound velocity along the 
chains. Such a length scale can be seen as the typical dis- 
tance above which the superfiuid phase fields of the two 
chains lock together. It is interesting to look at the ratio 
^a/d, where d is the mean "relevant distance" between 
particles. By relevant, we mean the distance between 
bosons: d ^ {2p)~^ for < p < 1/4 and between holes 
d ^ {1 — 2p)^^ for 1/4 < p < 1/2. There are two regimes 
of interest, (i) When ^ d, the system is effectively ID 
since the transverse phase fluctuations appear to decay 
over a distance much smaller than the inter-particle dis- 
tance d: in such a regime, we expect Kg ^ 1/2. (ii) In 
the opposite limit 3> d, the system will behave as two 
weakly interacting SF for which one expects Kg ^ 1. 

One can try to get a slightly more quantitative view 
on this, replacing K = 1 in Eq. (5.13), we get 



d 



.[(2p),(l-2p)] X 2sin(^p) 



t 



2/3 



(5.14) 
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FIG. 19: (Color online) Velocity Vc/2t, extracted from QMC 
simulations using Eq. (5.15), plotted against the density p for 
various values of the interchain hopping t±, and compared 
with two limiting cases (see text): Eq. (5.16) in the weak 
coupling limit (green full line), and Eq. (5.17) for the strong 
coupling limit (blue dashed line). 



Plotted in Fig. 18, this ratio displays an asymmetric 
behavior about p = 1/4 which is connected with the 
asymmetric behavior of Ks at small t±/t (Fig. 14). In 
particular, the much sharper vanishing of ^a/d close to 
half-filling than in the diluted limit explains the sharp 
decreasing to i^s = 1/2 at half- filling for small values 
of t±/t whereas at zero (or unit) filling Kg goes to 1/2 
much slower (see Fig. 14). One also sees that when t±/t 
increases, ^a/d becomes very small for all fillings. This 
clearly signals that the effective ID model obtained above 
in the other limit {t±/t ^ 1) becomes a good description 
and that the asymmetry tends to disappear, as already 
noticed. 



4- Critical velocity 

We finish this discussion with the analysis of the speed 
of sound V discussed above. This velocity can be in- 
terpreted as the SF critical velocity Vc associated to the 
linear branch in the excitation spectrum. Using the same 
hydrodynamic relations as above, we get for the velocity 



(5.15) 



QMC results are shown in Fig. 19 for Vc versus the density 
p in various limiting cases t± — OAt, t, 2t, IQt. For a 
very small interchain hopping we expect the decoupled 
chains result to hold for Vr, that is: 



Vc — 2t sin {np) if t± <^ t, 



(5.16) 



whereas in the very strong rung coupling limit, the effec- 
tive ID model Eq. (3.8) with correlated hopping terms 




density p 



FIG. 20: (Color online) Phase diagram for hard-core bosons 
on a two-leg ladder Eq. (1.1) in the plane density p - inter- 
chain hopping t±/{t + t±). The contour map for the Luttinger 
liquid parameter of the gapless symmetric mode Ks has been 
obtained from DMRG simulations. At p = 1/2 there is the 
gapped RMI (cyan line), and a SF phase for all other incom- 
mensurate fillings with Ks £ [1/2, 1]. 



yields a different velocity of excitations as follows: 



2i |sin(27r/3)| 



1 



4t 



(1 — p) cos(27rp) 



if t± > t. 



(5.17) 

The latter expression was obtained by following the 
calculation that led to the perturbative expression of Ks , 
equation 5.11. We computed the compressibility and 
superfluid stiffness of the effective model and used the 
hydrodynamic relation (5.15). The agreement between 
these two limiting cases and the QMC estimates for Vc 
are quite good, especially for the two limits t± = 0.4t 
and tj_ = lot. In between, for instance for t± = 2t, the 
shape of Vc{p) is clearly non- universal, and asymmetric. 



VI. DISCUSSIONS 

A. Phase diagram 

Before discussing further two important implications of 
our results, namely the effect of weak disorder on the SF 
phase and the consequences for experiments, we present 
the general phase diagram of the hard-core bosons model 
on a two-leg ladder Eq. (1.1). As visible in Fig. 20, in 
the plane [density p] — [t±/{t + t±)], a color map of the 
Luttinger parameter Kg of the gapless symmetric mode 
is shown. Except along the line p = 1/2, where there 
is the gapped RMI, the phase is SF with a continuously 
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varying Kg €E [1/2, 1]. The map has been obtained using 
DMRG results. All the limiting cases discussed above are 
visible here: Kg = 1/2 along p = 0, 1/2 and for all fillings 
when tj_/t — > oo; and Ks=l in the limit tj_/t — > 0. 



B. Disorder effects 

1. Weak disorder 

In this section we briefly discuss the pinning of a two- 
leg ladder boson gas on a random chemical potential. 
Following Ref.79, we consider a random chemical poten- 
tial along both chains. In the continuum, 



dx[Vi{x)pi{x) + V2{x)p2{x)]^ 



(6.1) 



We take Gaussian distributions for Vi , V2 and as- 
sume they are uncorrelated, that is Vi{x)Vi{x') = 
V2ix)V2{x') = DS{x - x') and Vi{x)V2ix') = 0. RG 
calculations lead to two very different results whether t± 
is very weak or very strong. In the case of uncoupled 
chains {t± — 0), the Bose gas is pinned on each chain 
by the 27rp Fourier component of the random potential 
as long as K < 3/2 and for very weak disorder one can 
estimate the localization length to be®°: 



ID 

loc 



2Da 



(6.2) 



The resulting localized phase is called a Bose glass 
(BG)^'*"'*^ which can be seen as an incoherent compress- 
ible ID Bose gas. Along the line a.t t±/t = 0, chains are 
decoupled and the BG phase consists of two independent 
strictly ID BG in which the density is pinned at 2ttp so 
that we use the denomination BG2-^p- Its localization 
length diverges aX K — 3/2, as quantum fluctuations de- 
stroy the pinning. 

On the contrary for very large transverse hopping 
induces long range phase coherence (in the antisymmetric 
mode), and the superfluid phase is stable against weak 
disorder provided Ks > 3/4. On our map of Ks (see 
Fig. 20), it appears that at moderate t±/t (roughly below 
1.5) there exists a region of densities for which Kg is 
indeed above the critical value of 3/4, making the system 
insensitive to a random chemical potential. In Fig. 21, 
we show the phase boundary Kg = 3/4 which separates 
regions where weak disorder D = 0^ is relevant (for Kg < 
3/4: BG4^p) or irrelevant (for Kg > 3/4: SF). In this 
limit, the transverse coherence length introduced in 
Sec. V C 2 is very short, and one expects a pinning of the 
Bose gas by the 47rp component of the random potential, 
hence the denomination BG^T^p. 



2. Finite disorder SF-BG transition 

One can make use of the criterion Kg > 3/4 for the sta- 
bility of the SF phase against weak disorder to plot three 
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FIG. 21: (Color online) Phase boundary at =3/4 (ob- 
tained from DMRG data) separating regions where weak dis- 
order is relevant (for Ks < 3/4: BG4,rp) or irrelevant (for 
Ks > 3/4: SF). Along the line at t±/t = 0, the chains are 
decoupled and two uncoupled BG at 2ttp are expected. 



representative cases at small, intermediate and strong 
transverse hopping. For instance at strong t±/t^ when Kg 
remains below the critical value of 3/4 for any filling, one 
cannot expect a stable SF phase at any finite disorder. 
This is what is represented in the right panel of Fig. 22 
for t±/t = 3. More precisely, we expect from the DMRG 
results for Kg that whenever t±/t > 1.5, since Kg < 3/4, 
weak disorder will immediately lead to localization. This 
result has to be contrasted with recent numerical stud- 
ies on disorder bosonic ladders®^ since we expect strong 
finite size effects (in particular at low fillings) to be an 
obstacle to a proper interpretation of numerical results. 
For intermediate transverse hopping t±/t ^ 1, we expect, 
as shown in the middle panel of Fig. 22 for t± = 1.2t, a 
stable SF phase for limited fillings but with a maximal 
critical disorder strength D^^'^^ not so small since one 
can naturally expect D"^^'^ ~ t^. For the third case at 
small t±/t (see Fig. 22 left), we expect a SF phase more 
extended in term of filling but with a reduced strength. 



with again D" 



t^. The asymmetry of Kg versus p 



will translate into a similar asymmetry for D^- At small 
densities, since Kg — >■ 1/2, we expect the BG phase to 
extend at weak disorder up to a finite critical concentra- 
tion pc- Close to the gapped RMI phase, while Kg 1/2 
much faster, we still expect a finite window (although 
much smaller) with pc < 1/2 where an intervening BG 
between RMI and SF phases is present. ^'^"^ Again this 
result has to be contrasted with the numerics obtained 
in Rcf.82 at small density and close to half-filling where 
finite size effects should be quite strong. 

In great contrast with the fermionic case where weak 
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FIG. 22: (Color online) Three typical situations for 2 coupled 
chains of hard-core bosons (small, intermediate and strong 
interchain hopping). On the top line, the Luttinger liquid pa- 
rameter of the gapless symmetric mode Ks is plotted (DMRG 
results) versus the density p. The critical line at Kl = 3/4 is 
also shown. Bottom: schematic phase diagram in the plane 
density p - disorder D displayed with two different phases, SF 
and BG at 47rp (see text). 

disorder would immediately localize the particles, hard- 
core bosons on a two-leg ladder provide a very simple 
model where a finite disorder phase transition between 
SF and BG phases is expected at incommensurate filling. 
This makes this model quite interesting in the view of 
large scale numerical studies to precisely investigate the 
SF-BG universality class, for which to the best of our 
knowledge, only a few studies exist. In particular, 
the question of the universal character of the Luttinger 
parameter right at the critical point for a finite disorder 
transition appears to be an important question for which 
two different scenarios have been put forward so far.^^'^^ 

C. Experimental consequences 

1. Coupled Base liquids in optical lattices 

A necessary requirement to probe our proposed phase 
diagram with cold atoms is first to obtain ladder-like 
structures from an optical lattice. This can be reached 
by realizing a lattice of double-wells which are achieved 
by superimposing two independent planar optical lattices 
with periodicity A and A/2. Such an array of double- 
wells has been proposed in Ref.86 and recently realized 
in Refs.21,87 with red-detuned lasers (see also Ref.88 for 
an extensive discussion on the realization of ladder-like 
potentials with lasers). A notable result of the present 
study is the phase diagram depicted in Fig. 2 where we 
predict the existence of a rung-Mott insulating phase at 
half filling whatever the transverse hopping parameter. 

However, we have shown that the charge gap scales 



as As ~ exp{—a.t/t±) which is exponentially small in the 
weak transverse hopping regime. Therefore for this phase 
to be observed, the condition ^ A^ must be satis- 
fied, where A^ ~ Hv/L and L is the longitudinal size 
of the cigar-shape trap and v the sound velocity along 
the chain. Since this size is finite, this puts a constraint 
on a minimum strength for the parameter t± in order to 
observe the rung-Mott insulator. Practically speaking, 
from Table I, a relatively large value oi t±/t ^ 1 is re- 
quired to make the Mott phase experimentally observable 
in cold atoms experiments. Such a Mott insulating phase 
is characterized by one atom per rung and should be de- 
tectable using single-atom resolved fiuorescence imaging 
as recently used for detecting a 2D Mott insulator. 

Technically, it seems rather straightforward to achieve 
the Tonks-Girardeau regime^^ for two coupled ID Bose 
liquids, as for instance in Ref.21. The system (of total 
length L) being trapped in a harmonic potential, there 
is a site-dependent chemical potential iJ,j = V{j — L/2)'^. 
Therefore, assuming a zero chemical potential at the trap 
center j — L/2, the Mott state having an energy gap A^ 
will extend over a finite length 




Using a realistic potential with V/t ~ 0.003,''^ the rung- 
Mott state will be detectable near the center of the trap 
over ^Mott ~ 20 sites for t±/t = 2. However, for smaller 
tj_/t the incompressible state will be practically unde- 
tectable since i?Mott decreases quite fast whereas in the 
same time, the correlation length ^ increases rapidly. 



2. Spin ladders in an external magnetic field 

Theoretically predicted some years ago,''^ field-induced 
Bose-Einstein condensation of triplet excitations in a 
magnetic insulator was first observed experimentally in 
the spin-1/2 dimer compound TlCuCla.^'^'^'* As reviewed 
in Ref.29, a quantitative theoretical description can be 
achieved through effective theories written in term of 
hard-core bosons. ^^'^^'^^ This is especially useful for sys- 
tems which, in zero field, have a singlet-triplet gap in the 
excitation spectrum, such as e.g. spin-1 chains, spin-1/2 
ladders, or more complex 3D arrays made of coupled spin- 
1/2 dimers for instance. While a true Bose condensed 
state is only possible for D> 1 systems, a spin ladder 
in an external field is nevertheless expected to achieve 
a Luttinger liquid, with SF quasi-long-range order. Re- 
cently, an excellent experimental realization for such a 
Luttinger liquid physics was discovered in the ladder ma- 
terial (C5Hi2N)2CuBr4 (BPCB)27,28,96,97 foj. ^-^^^-^ ^ 

rect computation of the Luttinger parameters K and Vc 
has been performed'^ ^ using various numerical and ana- 
lytical techniques. BPCB is well described by a spin-i 
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antiferromagnetic Heisenberg model on a ladder 

n = J J2 • S,,,+i + ^ Si,, • S2,j (6.4) 

with 7 = J/ J±. Nuclear magnetic resonance (NMR) and 
inelastic neutron scattering (INS) both give 7 ~ 0.25,^^'^^ 
which puts this material in the class of strong rung cou- 
plings. As performed in Sec. Ill A 2 and in the work by 
Bouillot et al.,^^ one can write down an effective model 
with hard-core bosons on a chain with correlated hopping 
and nearest-neighbor interactions. Here, interactions are 
repulsive because of the antiferromagnetic exchange that 
is only slightly shielded by the second-order attractive 
interactions. Although the perturbative approach devel- 
opped in Sec. V C 2 is unable to describe quantitatively 
the evolution of K with the magnetization per rung - 
mainly because of the large value of the repulsion in the 
effective model - it reproduces some of its features. It 
correctly gets the limits X — 1 as and 1. It also 

sheds light on the competing effects of correlated hop- 
pings and repulsive nearest-neighbor interactions. Cor- 
related hopping tends to renormalize K upwards, which 
explains the departure from the values of K obtained 
from the Bethe ansatz solution of the XX Z chain. Its 
effects are more pronounced for < 1/2, since next- 
nearest neighbor hopping is shielded as the filling of the 
triplon band approaches 1. The case of BPCB is interest- 
ing in the sense that it lies on the strong rung coupling 
side with 7 ~ 1/4 where the Luttinger liquid physics in 
the strong coupling limit can be tested. In the other 
limit (7 ^ 1) a Luttinger liquid regime is also expected, 
as recently observed in the compound DIMPY by Hong 
and co-workers'^'^ with 7^2. Nevertheless, as discussed 
previously for the case of hard-core bosons, we expect 
a strong anisotropy in the Luttinger parameters, as ob- 
served for the velocity in DIMPY. 

One can also mention two other interesting materials: 
IPA-CuCla and BiCu2P06 which both have a sizeable 
spin gap As - 10 T for the former^^.^s and ~ 22 T 
for the latter. They both display interesting ladder- 
like features with some remarkable peculiarities, namely 
a strong ferromagnetic diagonal coupling in the case 
of IPA-CuCla^"^ and a presumably non-negligible intra- 
chain frustration in the case of BiCu2P06. These two ef- 
fects may lead to a non-trivial field-dependent Luttinger 
liquid physics. Moreover, the exploration of disorder ef- 
fects have recently started in these materials^^^"^"'^ where 
one expects a Bose glass physics to occur when the gap is 
closed by the external magnetic field. An interesting sit- 
uation would be to observe a SF-BG transition induced 
by the external field. Of course a dimensional crossover 
is also expected at low temperature. 

VII. CONCLUSIONS 

In this paper, we have presented an exhaustive study 
of the behavior of hard-core bosons on a two-leg lad- 



der. The two-leg ladder is a simple minimal geome- 
try that allows for particle exchange, in contrast with 
strictly one dimensional chains (limited to nearest neigh- 
bor hopping). Therefore, particle statistics should affect 
the physical behavior of the system. Using a combination 
of several analytical and numerical techniques, we have 
indeed shown explictly that the phase diagram of hard- 
core bosons on a two-leg ladder differs significantly from 
the one of free fermions as summarized in Fig. 2. The 
most stringent difference occurs at half-filling and small 
transverse hopping where hard-core bosons enter a rung- 
Mott insulating phase while free fermions are still in a 
metallic phase. Outside half-filling, hard-core bosons are 
in a superfluid phase which is characterized by the Lut- 
tinger liquid exponent Kg which varies considerably with 
the density between 1/2 and 1. This result is well summa- 
rized in Fig. 20 where Kg is plotted as a function of the 
density and the transverse hopping amplitude. The Lut- 
tinger liquid parameter governs all correlation functions 
but also tells us whether the system might be sensitive 
to weak disorder or not. Our quantitative determina- 
tion of Kg therefore allowed us to draw qualitatively a 
phase diagram for hard-core bosons on a two-leg ladder 
with weak disorder (see Figs 21 and 22). In great con- 
trast with the fermionic case where weak disorder would 
immediately localize the particles, hardcore bosons on a 
two-leg ladder provide a very simple model where a fi- 
nite disorder phase transition between a superfluid and 
Bose glass phases is expected at incommensurate filling. 
This makes this model quite interesting in this respect. 
A more elaborated analysis of the effects of disorder on 
the behavior of hard-core bosons on a two-leg ladder will 
be the subject of a subsequent analysis. 
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Appendix A: Gap at half-filling 

We outline here the perturbative calculation to second 
order in t/t± that led to the expression of the size of 
the plateau at half-filling. We start from the limit of 
decoupled rungs i/t_L =0. The ground state Hamiltonian 
reads 

^o = -t±E(^tA,+h.c.). (Al) 
j 

Four states are available on a given rung j: an 
empty state \0)j, two 1-particle state |l±)j = 
(|j, 1) ± |j,2)) /V2, and a 2-particle state \2)j. At half- 
filling, N = L, the ground state is a tensor product of 
L symmetric states, = |l-l-, 1-I-, . . . , 1-I-) and the 
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ground state energy is Ef'^ = -t^L. Now we turn on 
longitudinal hopping, with i <C i_L. The perturbation 
Hamiltonian is 

Hi = -i^(5l_^.5,_^.+i+h.c.). (A2) 

Acting on li*-*^-*), Hi creates 2L particle-hole excitations 
on nearest-neighbor rungs, 



state). This formula has a very intuitive form: Hi cre- 
ates particle-hole excitations, that are selected by Qo- 
Then Hi acts a second time on this excited state. A 
non zero correction is obtained only if this second pro- 
cess takes the system back to the ground-state subspace. 
Here the excited states have an energy -~t±{L — 2) and 

Qo[E^L^ - Ho]-^Qo can be replaced by Qo/{2tj_). The 
second order correction reads: isj^^'' = —2Lt^/{2t±), and 
we find for the energy of the L-particle state. 



Hi\L(^^) = -tY: 



2 0,...,1+) 

0) 



.,0 2,...,: 

u) 



(A3) 



There is no first-order correction to the ground-state en- 
ergy since (L(o)|iJi|L(o)) = 0. The second-order correc- 
tion is of the form^^ 



E 



(2) 



{L^'^\HiQo^ goifi|i'°^), (A4) 



4"^ - ^0 



where Qq is a projector, that projects out of the ground- 
state subspace (composed here of one non-degenerate 



Er 



-tj_L - 2L 



2t\ 



(A5) 



We now turn to the states with L + 1 particles. They are 
L times degenerate and are best written in a momentum 
representation. 



(A6) 

Their ground-state energy is = —t±{L—l). Hi can 

either move the extra particle around (to next-nearest 
rungs) or create 2{L — 1) particle-hole excitations: 



iJi|(L + l)i°') = 



t 

/2L 



^e"^-'- ( |l+,l-f,...,l+,, 2 ,,1- 



(j+i) 



E 



,1- 



u-i) 



2,0 ,21+,. 

(i-2) 0-1) (j) 



1+) 



1+, . . . , 2 , 0, . . . , 2 1+, . . . , 1+) + 1+, . . . , 2, . . . , 2 1+, . . . , 1+) 

(p) u) (p) (i) 



U) (j + l) (j+2) 



(A7) 



The first line is the result of moving the extra-particle to 
neighboring rungs. It brings a first-order correction to 
the {L + l)-particle state energy, £;W^(fc) = -2icos(fc). 
The second and third lines involve excited states that will 
bring a correction to second order: 

E^^l, = {{L + l)(°)|HiQ^_J Q^Hi\{L + 

^L+i - Ho 

Excited states have here an energy of —t±{L — 3) and we 
can again replace Qo[E^^^^ — Hq]~^Qo by Qo/i'2t±)- Act- 
ing a second time Hi will shift the energy by a constant 
(acting on the states on the second line) as it did for the 
energy of the L-particle state. It will also move the extra 
particle to next-nearest rungs (acting on the states on 
the third line). The second order correction finally reads 
E'-j^l^ik) = -2{L - 2)t^/{2t±) - 2tV(2tj_)cos(2fc), and 



the corrected energy for the L -\- 1-particle state is: 
EL+iik) = -tj_{L-l)-2tcos{k) 

- 2(i-2)^-2|^cos(2/c). (A8) 

Note that we could do the exact same work to compute 
the energy of the state with L ~ 1 particles and would 
find the same energy. Thus, the size of the plateau at 
half-fining is given by = 2[El+i{0) - El]: 

2f 

As = 2i_L-4tH . (A9) 

tj_ 

Appendix B: Effective model 

One can use a somewhat different method to obtain 
the effective Hamiltonian to second order in perturbation 
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theory for incommensurate fillings. ^^'^^ We start from the 
original model that we cast into, H = Hq + "Hi, with: 

-Hi = -tY^ (^Ij^lj+i + h-c.) - (m - A*c) ^t^d 
j.e 3,1 

"Ho = -t± Y (^1 j4j + ^■'^■) ^ X! "^J 

(Bl) 

and /ic = i_L. For the ground state Hamiltonian 'Ho, 
states and \2)j have the same energy = —2t±. 

Let's call Pq the projector on the fundamental subspace, 
made of 2^ states. In this subspace and \2)j are 

the only states allowed on a given rung j. Let us call 
Qo the projector on the complementary subspace. The 
effective Hamiltonian to second order is: 

Heff = PoHiPo + PoniQo[E - Ho]-^QoniPo, (B2) 

with E the eigenvalue of degenerate subspace under con- 
sideration. Here, given the degenerate ground state and 
the form of Hi, the virtual excited state will be a state 
with an empty rung. Then, in everything that follows, 
Qo[E ~ Hor'Qo - Qo[-2t^L + 2t^{L - l)]-'Qo^ - 
--Qo/{2t±). In the first order contribution the action 
of the projectors Pq can be enforced by using spin-1/2 
operators instead of bosons: 



^3,1 = bj,2 



72^^ 
72^^ 



1 














+ 1 


2 







(B3) 
(B4) 
(B5) 



Then, 



E 

3 



3 

■ + 2 
(B6) 

Now we need to compute PoHiQaHiPo- We have seen 
above in the previous calculation that the only way to 
create an excited state is to let V-i empty a singly oc- 
cupied rung. Let j be such a rung. Then the term of 
interest in QqHiPq is of the form 



(1 - Tlj- 1 - nj- 2)(6]+i,i^- 1 + &]+i,2^,2)-Pc 



(B7) 



Now one must act with PqHiQq on the latter term and 
bring back the system into the original subspace. As we 
have already seen, there are two ways of doing so: 



-^l2^i+i,2)K+ia+"i+i,2- 



1) 



a) Act with ^o(6j,ifcj+i,i 

and get back the initial state. The corresponding con- 
tribution to the second order correction is an attractive 
nearest-neighbor interaction of the form: 

_ ^ 
t± 



E 



(B8) 



where we have used the spin 1/2 operators to enforce 
the projection. 

b) Act withPo(fe],i&-,_i,i+^],2^i-i,2)("-i-ia+"-i-i,2-l)- 
This generate a next nearest neighbor hopping term of 
the form: 



2tj^ 



E4+1 



'3-1- 



(B9) 



A Jordan- Wigner transformation transforms the latter 
Hamiltonian into one of a band of spinless fermions with 
next-nearest neighbor hopping and attractive interac- 
tions: 



H„(f = —t 



t± 



^ [c]cj+i -t- H.c. -(pt-//J^(nj + l) 

3 j 

3 

^(l-nj)(l-nj+i) 



(BIO) 



Appendix C: Calculation of Ks and Vc in the strong 
coupling regime 

We briefly outline here the calculation that led to the 
perturbative expressions (5.11) and (5.17) for Kg and 
Vc- We start from the effective Hamiltonian of equation 
(3.8), derived in section III A 2, and extract its Luttinger 
liquid parameters v and K. The low-energy form of the 
Hamiltonian is^°^: 



H 



h 

2^ 



dx[vji\7ef+VN{V^f], (CI) 



where the density and phase stiffness vn and vj are re- 
lated to the ground-state energy as^^: 



VN 



L_ f d'Ecs ] 
Ln r d^EGs \ 



(C2) 
(C3) 



Here, 95 is a twist angle imposed on the boundary con- 
ditions as c|^^ = e*'^c] and iVp is the number of par- 
ticles. Following Cazalilla,^^ the ground-state energy is 
computed perturbatively to first order in t/t^_ as Eqs = 
(Hq)o + (Hi)o, with 



Hn= -t 



h.c 



H, = - 



(C4) 



2t7 



3 



(C5) 
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Using Wick's theorem we find: 

Egs= -tL{h+f^i) 



2t± 



L [(1 - /o)(/2 + f-2) + (/i - I-if + 2fl] 



(C6) 

with /p = (c]^pCj)o- Note that fo = p = Nq/L is the 
density of spinless fermions. Leaving aside the algebra 
we find, when taking L — >■ oo: 

2t2 2^2 
vn = 2t sin(7rp) H cos(27r/5) 



2t2 

+ ^(1-/5) sin(2^p), 
iij = 2ts'm(np) H sin (tt/j) H 



(C7) 



(1 — p) sin(27rp). 

(C8) 



Finally, we identify the Luttinger parameters v and K 
by using the relations vK — vj and u/K — vm (note the 
connexion with equations (5.3) and (5.15) in the main 
text) and find: 



K = 1 
V = 2tsai{'Kp) 



2t sin(7rp) 



2t 

1 H {1 — p) cos(7rp) 



(C9) 
(CIO) 



We argue in the main text that for large values of t^/t 
the correct identification for the original Luttinger pa- 
rameters is Kg = K /2 and Vc = v. 



XXZ) for which one can compute exactly the SF density 
using the Bethe Ansatz.^^'^^ This model, governed by the 
Hamiltonian 



n 



t-v 



= t E + h-c-] + E Vn.^^+i , (Dl) 



is very well known. In particular, there is a SF- 
Insulator transition where, in the thermodynamic limit, 
the SF density jumps between 1/4 and 0^^'''^ at the 
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FIG. 23: (Color online) Betlie Ansatz results for the SF den- 
sity of the integrable t — V chain in the repulsive regime 
obtained for lengths L — 8, ■ ■ ■ , 1024 and various repulsion 
strengths V. Once the data collapse is obtained, the result- 
ing correlation length ^ is displayed in the inset (red squares) 
and successfully compared to the exact results from Baxter^ 
(green dashed curve). 



Appendix D: Correlation length from data collapse 
in the t — V model 

In this appendix, we show the usefulness of the "scal- 
ing plot" technique to obtain the correlation length. As 
already used in a context of disordered chains to de- 
termine the localization length, it turns out to be a quite 
efficient method to access (i) the value of the critical cou- 
pling, and (ii) the way the correlation (or localization) 
length diverges at criticality, even when the numerically 
accessible system sizes are much smaller than the ac- 
tual correlation length. To illustrate this technique on a 
well-controlled example, we briefly look at the hard-core 
bosonic t—V model (or equivalently fermionic or spin-1/2 



Kosterlitz-Thouless transition point V = 2t. On the in- 
sulating side V > 2t, the correlation length ^, computed 
exactly by Baxter, diverges exponentially. For finite 
length chains, the SF density vanishes ~ exp(— i/f) 
while it is finite on the other side of the transition. Us- 
ing the exact Bethe Ansatz solution of this model with 
twisted boundary conditions, the SF density is com- 
puted^^ and the scaling plot analysis performed above for 
the ladder system is repeated in Fig. 23 where one can 
see the two regimes of Eq. (4.12). The correlation length 
extracted from the collapse is in very good agreement 
with the exact result^^^ as shown in the inset of Fig. 23. 
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